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O ■ I. INTRODUCTION 

o : 
o . 

CN ■ A. On approximation methods in general relativity 

5h ! 
Ph. 

^ . Let us declare that the most important devoir of any physical theory is to draw firm 

^ . predictions for the outcome of laboratory experiments and astronomical observations. Un- 
fortunately, the devoir is quite difficult to fulfill in the case of general relativity, essentially 
^ . because of the complexity of the Einstein field equations, to which only few exact solutions 
^ ! are known. For instance, it is impossible to settle the exact prediction of this theory when 
O ! there are no symmetry in the problem (as is the case in the problem of the gravitational 
^ ! dynamics of separated bodies). Therefore, one is often obliged, in general relativity, to resort 
O \ to approximation methods. 

^ ! It is beyond question that approximation methods do work in general relativity. Some of 

qh| the great successes of this theory were in fact obtained using approximation methods. We 
^ ! have particularly in mind the test by Taylor and collaborators f^-^H regarding the orbital 
?^^! decay of the binary pulsar PSR 1913+16, which is in agreement to within 0.35% with the 
• w^ ', general- relativistic post- Newtonian prediction. However, a generic problem with approxi- 
^ \ mation methods (especially in general relativity) is that it is non trivial to define a clear 
5^ \ framework within which the approximation method is mathematically well-defined, and such 
that the results of successive approximations could be considered as theorems following some 
precise (physical and/or technical) assumptions. Even more difficult is the problem of the 
relation between the approximation method and the exact theory. In this context one can 
ask: What is the mathematical nature of the approximation series (convergent, asymptotic, 
. . .)? What its "reliability" is (i.e., does the approximation series come from the Taylor ex- 
pansion of a family of exact solutions)? Does the approximate solution satisfy some "exact" 
boundary conditions (for instance the no- incoming radiation condition)? 

Since the problem of theoretical prediction in general relativity is complex, let us distin- 
guish several approaches (and ways of thinking) to it, and illustrate them with the example 
of the prediction for the binary pulsar. First we may consider what could be called the "phys- 
ical" approach, in which one analyses the relative importance of each physical phenomena at 
work by using crude numerical estimates, and where one uses only the lowest-order approxi- 
mation, relating if necessary the local physical quantities to observables by means of balance 
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equations (perhaps not well defined in terms of basic theoretical concepts). The physical 
approach to the problem of the binary pulsar is well illustrated by Thorne in his beautiful 
Les Houches review [Q (see also the round table discussion moderated by Ashtekar 0): one 
derives the loss of energy by gravitational radiation from the (Newtonian) quadrupole for- 
mula applied formally to point-particles, assumed to be test-masses though they are really 
self-gravitating, and one argues "physically" that the effect comes from the variation of the 
Newtonian binding energy in the center-of-mass frame - indeed, on physical grounds, what 
else could this be (since we expect the rest masses won't vary)? The physical approach 
yields the correct result for the rate of decrease of the period of the binary pulsar. Of 
course, thinking physically is extremely useful, and indispensable in a preliminary stage, 
but certainly it should be completed by a solid study of the connection to the mathematical 
structure of the theory. Such a study would a posteriori demote the physical approach to 
the status of "heuristic" approach. On the other hand, the physical approach may fall short 
in some situations requiring a sophisticated mathematical modelling (like in the problem 
of the dynamics of singularities), where one is often obliged to follow one's mathematical 
rather than physical insight. 

A second approach, that we shall qualify as "rigorous", has been advocated mainly by 
Jiirgen Ehlers (see, e.g., 0). It consists of looking for a high level of mathematical rigor, 
within the exact theory if possible, and otherwise using an approximation scheme that we 
shall be able to relate to the exact theory. This does not mean that we will be so much 
wrapped up by mathematical rigor as to forget about physics. Simply, in the rigorous 
approach, the prediction for the outcome of an experiment should follow mathematically 
from first theoretical principles. Clearly this approach is the one we should ideally adhere 
to. As an example, within the rigorous approach, one was not permitted, by the end of 
the seventies, to apply the standard quadrupole formula to the binary pulsar. Indeed, as 
pointed out by Ehlers ei a/ [0, it was not clear that gravitational radiation reaction on a 
self-gravitating system implies the standard quadrupole formula for the energy flux, notably 
because computing the radiation reaction demands a priori three non-linear iterations of the 
field equations P], which were not fully available at that time. Ehlers and collaborators 
remarked also that the exact results concerning the structure of the field at infinity (notably 
the asymptotic shear of null geodesies whose variation determines the flux of radiation) were 
not connected to the actual dynamics of the binary. 

Maybe the most notable result of the rigorous approach concerns the relation between 
the exact theory and the approximation methods. In the case of the post-Newtonian ap- 
proximation (limit c —>■ oo), Jiirgen Ehlers has provided with his frame theory ||9|-[TT| a 
conceptual framework in which the post-Newtonian approximation can be clearly formu- 
lated (among other purposes). This theory unifies the theories of Newton and Einstein into 
a single generally covariant theory, with a parameter 1/c taking the value zero in the case of 
Newton and being the inverse of the speed of light in the case of Einstein. Within the frame 
theory not only does one understand the limit relation of Einstein's theory to Newton's, 
but one explains why it is legitimate when describing the predictions of general relativity 
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to use the common-sense language of Newton (for instance thinking that the trajectories of 
particles in an appropriately defined coordinate system take place in some Euclidean space, 
and viewing the coordinate velocities as being defined with respect to absolute time). It was 
shown by Lottermoser |]r2| that the constraint equations of the (Hamiltonian formulation of 
the) Ehlers frame theory admit solutions with a well-defined post-Newtonian limit. Further 
in the spirit of the rigorous approach, we quote the work of Rendall on the definition 
of the post-Newtonian approximation, and the link to the post-Newtonian equations used 
in practical computations. (See also [l^Jl^ for an attempt at showing, using restrictive 
assumptions, that the post-Newtonian series is asymptotic.) 

The important remarks of Jiirgen Ehlers et al on the applicability of the quadrupole 
formula to the binary pulsar stimulated research to settle down this question with (al least) 
acceptable mathematical rigor. The question was finally answered positively by Damour 
and collaborators who obtained in algebraically closed form the general- relativistic 

equations of motion of two compact objects, up to the requisite 5/2 post-Newtonian order 
(2.5PN order or 1/c^) where the gravitational radiation reaction force appears. This ex- 
tended to 2.5PN order the work at IPN of Lorentz and Droste [^, and Einstein, Infeld 
and Hoffmann The net result is that the dynamics of the binary pulsar as predicted 
by (post-Newtonian) general relativity is in full agreement both with the prediction of the 
quadrupole formula, as derived earlier within the "physical" approach, and with the obser- 
vations by Taylor et al (see [22| for discussion). 

Motivated by the success of the theoretical prediction in the case of the binary pulsar 
T^-^^ , we shall try to follow in this article the spirit of the "rigorous" approach of Jiirgen 



Ehlers, notably in the way it emphasizes the mathematical proof, but we shall also differ 
from it by a systematic use of approximation methods. This slightly different approach rec- 
ognizes from the start that in certain difficult problems, it is impossible to derive a physical 
result all the way through the exact theory without any gap, so that one must proceed 
with approximations. But, in this approach, one implements a mathematically well-defined 
framework for the approximation method, and within this framework one proves theorems 
that (ideally) guarantee the correctness of the theoretical prediction to be compared with 
experiments. Because the comparison with experiments is the only thing which matters in 
fine for a pragmatist, we qualify this third approach as "pragmatic". 

In this article we describe the pragmatic approach to the problem of gravitational radia- 
tion emitted by a general isolated source, based on the rigorous post-Minkowskian iteration 
of the field outside the source |23|, and on the general connection of the exterior field to 
the field inside a slowly-moving source ||2^j25[] . Note that for this particular problem the 
pragmatic approach is akin to the rigorous one in that it permits to establish some results 
on the connection between approximate and exact methods. For instance it was proved by 
Damour and Schmidt p6[ (see also [^,^]) that the post-Minkowskian algorithm generates 
an asymptotic approximation to exact solutions, and it was shown [|29| that the solution 
satisfies to any order in the post-Minkowskian expansion a rigorous definition of asymptotic 
flatness at future null infinity. However it remains a challenge to analyse in the manner of 
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the rigorous approach the relation to exact theory of the whole formalism of p3|-p5|,p£ 



By combining the latter post-Minkowskian approximation and a post-Newtonian expan- 
sion inside the system, it was proved (within this framework of approximations) that the 
quadrupole formula for slowly-moving, weakly-stressed and self-gravitating systems is cor- 
rect, even including post-Newtonian corrections [BDl; and idem for the radiation reaction 



forces acting locally inside the system, and for the associated balance equations ||31| , |32 



These results answered positively Ehlers' remarks in the case of slowly-moving extended 
(fluid) systems. However we are also interested in this article to the application to binary 
systems of compact objects modelled by point-masses. Indeed the latter sources of radiation 
are likely to be detected by future gravitational-wave experiments, and thus concern the 
pragmatist. We shall see how one can address the problem in this case. (When special- 
ized to point-mass binaries, the results on radiation reaction pT| , |32| are in agreement with 
separate work of Iyer and Will | 5B|J5^ .) For other articles on the problem of gravitational 
radiation from general and binary point-mass sources, see p5H39| . 



B. Field equations and the no-incoming radiation condition 

The problem is to find the solutions, in the form of analytic approximations, of the 
Einstein field equations in iR^, 

and thus also of their consequence, the equations of motion of the matter source, V^T^'^ = 0. 
Throughout this work we assume the existence and unicity of a global harmonic (or de 
Bonder) coordinate system. This means that we can choose the gauge condition 

duh^"" = ; h^"^ = y/^g^" - T]^" , (1.2) 

where g and g^'^ denote the determinant and inverse of the covariant metric g^,^, and where 
1]^'^ is an auxiliary fiat metric [i.e. rj^'^ = diag(— 1, 1, 1, 1) = rj^^,]. The Einstein field equations 
(1) can then be replaced by the relaxed equations 

nh>"' = -^r'^" , (1.3) 

where the box operator is the flat d'Alembertian, □ = = rj'^'^d^du, and where the source 
term is the sum of a matter part and a gravitational part, 

^ I l^M^- + Af^:^ . (1.4) 

In harmonic coordinates the field equations take the form of simple wave equations, but 
whose source term is actually a complicated functional of the gravitational field h^'^; notably 
the gravitational part depends on h^'^ and its first and second space-time derivatives: 
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- g^'gardxh'^^dph'^' - g'^'^g^rdxh^^dph'^^ + gp^g^^d^h^^Pdrh'''' 

+ ^{2g^'g'"' - g^'gni'^gxrge. - gregx.)dph^^d,K'' . (1.5) 

The point is that A'*'^ is at least quadratic in /i, so the relaxed field equations (3) are very 
naturally amenable to a perturbative non-linear expansion. As an immediate consequence 
of the gauge condition (2), the right side of the relaxed equations is conserved in the usual 
sense, and this is equivalent to the equations of motion of matter: 

d^T^" = ^ V^T'^'^ = . (1.6) 

We refer to r'''^ as the total stress-energy pseudo-tensor of the matter and gravitational fields 
in harmonic coordinates. Since the harmonic coordinate condition is Lorentz covariant, r'*'^ 
is a tensor with respect to Lorentz transformations (but of course not with respect to general 
diffeomorphisms) . 

In order to select the physically sensible solution of the field equations in the case of 
a bounded system, one must choose some boundary conditions at infinity, i.e. the famous 
no-incoming radiation condition, which ensures that the system is truly isolated (no radi- 
ating sources located at infinity). In principle the no- incoming radiation condition is to be 
formulated at past null infinity J~ . Here, we shall simplify the formulation by taking advan- 
tage of the presence of the Minkowski background 77^,^ to define the no-incoming radiation 
condition with respect to the Minkowskian past null infinity J-^. Of course, this does not 
make sense in the exact theory where only exists the metric g^^^ and where the metric 77^,^ 
is fictituous, but within approximate (post-Minkowskian) methods it is legitimate to view 
the gravitational field as propagating on the fiat background r]ni,, since -q^i, does exist at any 
finite order of approximation. 

We formulate the no-incoming radiation condition in such a way that it suppresses any 
homogeneous, regular in IR^, solution of the d'Alembertian equation Uh — 0. We have at 
our disposal the Kirchhoff formula which expresses /i(x', t') in terms of values of /i(x, t) and 
its derivatives on a sphere centered on x' with radius p = |x' — x| and at retarded time 
t = t' - p/c: 
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d , Id 



dp c dt 



where dQ is the solid angle spanned by the unit direction (x — x') / p. From the Kirchhoff 
formula we obtain the no-incoming radiation condition as a limit at J'^, that is r — > -|-oo 
with t + r/c =const (where r = |x|). In fact we obtain two conditions: the main one. 



lim 

r— » + oo 
t+r /c=const 



r\ 1 

^irhn + -^Arhn (x,i) = 0, (1.8) 
or cot 
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and an auxiliary condition, that rd\h^" should be bounded at coming from the fact 
that p in the Kirchhoff formula (7) differs from r [we have p = r — x'.n + 0(l/r) where 
n = x/r]. 

In fact, we adopt in this article a much more restrictive condition of no- incoming radia- 
tion, namely that the field is stationary before some finite instant — T in the past: 

t<-T^ |[/^^'^(x,t)]=0. (1.9) 

In addition we assume that before — T the field h^'^{x.) is of order 0{l/r) when r +00. 
These restrictive conditions are imposed for technical reasons following , since they allow 
constructing rigorously (and proving theorems about) the metric outside some time-like 
world tube r = |x| > 7^. We shall assume that the region r > TZ represents the exterior of 
an actual compact-support system with constant radius d <7l [i.e. d is the maximal radius 
of the adherence of the compact support of T'^'^(x, t), for any time t]. 

Now if h'^'^ satisfies for instance (9), so does the pseudo-tensor r^*^ built on it, and then 
it is clear that the retarded integral of t^'^ satisfies itself the same condition. Therefore one 
infers that the unique solution of the Einstein equation (3) satisfying the condition (9) is 

= ^^D-V/^- , (1.10) 

where the retarded integral takes the standard form 

(□^V)(x,t)^-i- I ^^r(x',t-|x-x'|/c) . (1.11) 

Notice that since t'^'^ depends on h and its derivatives, the equation (10) is to be viewed 
rather as an integro-differential equation equivalent to the Einstein equation (3) with no- 
incoming radiation. 



C. Method and general physical picture 

We want to describe an isolated system, for instance a "two-body system" , in Einstein's 
theory. We expect (though this is not proved) that initial data sets g^u, dtg^u, p, v satisfying 
the constraint equations on the space-like hypersurface t = exist, and that this determines 
a unique solution of the field equations for any time t, which approaches in the case of 
two bodies a "scattering state" when t —00, in which the bodies move on unbound 
(hyperbolic-like) orbits. We assume that the space-times generated by such data admit a 
past null infinity JT" (or, if one uses approximate methods, J'^) with no incoming radiation. 
(Note that in a situation with initial scattering the field might not satisfy the rigorous 
definitions of asymptotic flatness at JT""; see [iD-Q.) The point to make is that in this class 



of space-times there is no degree of freedom for the gravitational field (we could consider 
other situations where the motion is influenced by incoming radiation). 
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Both our technical assumptions of compact support for the matter source (with constant 
radius d) and stationarity before the time — T contradict our expectation that a two-body 
system follows an unbound orbit in the remote past. We do not solve this conflict but 
argue as follows: (i) these technical assumptions permit to derive rigorously some results, 
for instance the expression [given by (52) with (56) below] of the far-field of an isolated past- 
stationary system; (ii) it is clear that these results do not depend on the constant radius 
d, and furthermore we check that they admit in the "scattering" situation a well-defined 
limit when — T — oo ; (iii) this makes us confident that the results are actually valid for 
a more realistic class of physical systems which become unbound in the past and are never 
stationary (and, even, one can give a posteriori conditions under which the limit — T —oo 
exists for a general system at some order of approximation). 

Suppose that the system is "slowly-moving" [in the sense of (12) below], so that we can 
compute the field inside its compact support by means of a post-Newtonian method, say 
h^^ = h^'^ where the overbar refers to the /orma/ post-Newtonian series. The post-Newtonian 
iteration (say, for hydrodynamics) is not yet defined to all orders in 1/c, but many terms 



are known: see the works of Lorentz and Droste 20 , Einstein, Infeld and Hoffmann 21 



Fock [44], Chandrasekhar and collaborators [45-47], Ehlers and followers [48-54], and many 
other authors []55|- p^ , p0| ,^ . 

On the other hand, outside the isolated system, the field is weak everywhere and it 
satisfies the vacuum equations. Therefore, the equations can be solved conjointly by means 
of a weak-field or post-Minkowskian expansion {G — >■ 0), and, for each coefficient of 
in the latter expansion, by means of a multipole expansion (valid because we are outside). 
The general Multipolar-post-Minkowskian (MPM) metric was constructed in p3| , pU[ as a 
functional of two sets of "multipole moments" Mi(t) and Siit) which were left arbitrary at 
this stage (i.e. not connected to the source). The idea of combining the post-Minkowskian 
and multipole expansions comes from the works of Bonnor [|5^ and Thorne [|60|. We denote 
by /igxt = M.{h^'^) the exterior solution, where M. stands for the multipole expansion (as it 
will turn out, the post-Minkowskian expansion appears in this formalism to be somewhat 
less fundamental than the multipole expansion). 

The key assumption is that the two expansions h!;^ = h^" and hZt = -M(/i^^) should 
match in a region of common validity for both the post-Newtonian and multipole expansions. 
Here is where our physical restriction to slow motion plays a crucial role, because such an 
overlap region exists (this is the so-called exterior near-zone) if and only if the system 
is slowly-moving. The matching is a variant of the well-known method of matching of 
asymptotic expansions, very useful in gravitational radiation theory []6Tl-|65|, ^ , |6^ ,p^ . It 
consists of decomposing the inner solution into multipole moments (valid in the outside), 
re-expanding the exterior solution in the near zone (r/c 0), and equating term by term 
the two resulting expansion series. From the requirement of matching we obtain in [^ , and 
review in Sections 2 and 3 below, the general formula for the multipole expansion }A{h^^) 
in terms of the "source" multipole moments (notably a mass-type moment II and a current- 
type Jl), given as functionals of the post- Newtonian expansion of the pseudo-tensor, i.e. 
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r^*^. [The previous moments Ml and Sl (referred below to as "canonical") are deduced from 
the source moments after a suitable coordinate transformation.] In addition the matching 
equation determines the radiation reaction contributions in the inner post-Newtonian metric 

To obtain the source multipole moments in terms of basic source parameters (mass 
density, pressure), it remains to replace t^^ by the result of an explicit post-Newtonian 
iteration of the inner field. This was done to IPN order in |]30|j66[| , then to 2PN order in 



2^ , and the general formulas obtained in permit recovering these results. See Section 



6. On the other hand, if one needs the equations of motion of the source, simply one inserts 
the post-Newtonian metric into the conservation law dpT^^ = 0. (Note that we are speaking 
of the equations of motion, which take for instance the form of Euler-type equations with 
many relativistic corrections, but not of the solutions of these equations, which are typically 
impossible to obtain analytically.) 



From the harmonic coordinates, one can perform to all post-Minkowskian orders ||29 



a coordinate transformation to some radiative coordinates such that the metric admits a 
far-field expansion in powers of the inverse of the distance R (without the powers of In R 
which plague the harmonic coordinates). Considering the leading order 1/R one compares 
the exterior metric, which is parametrized by the source moments (connected to the source 
via the matching equation), to the metric defined with "radiative" multipole moments, say 
Ul and V^. This gives Ul and Vl in terms of the source moments, notably II and Jl, and a 
fortiori of the source parameters. This solves, within approximate methods, the problem of 
the relation between the far field and the source. The radiative moments have been obtained 
with increasing precision reaching now 3PN pU|-[7I[|, as reviewed in Section 5. 

The previous scheme is developed for a general description of matter, however restricted 
to be smooth (we have in mind a general "hydrodynamical" T^"). Thus the scheme a priori 
excludes the presence of singularities (no "point-particles" or black holes), but this is a 
serious limitation regarding the application to compact objects like neutron stars, which can 
adequately be approximated by point-masses when studying their dynamics. Fortunately, 
the formalism is applicable to a singular T'^^ involving Dirac measures, at the price of a 
further ansatz, that the infinite self-field of point-masses can be regularized in a certain 
way. By implementing consistently the regularization we obtain the multipole moments and 
the radiation field of a system of two point-masses at 2.5PN order D72|[75[ , as well as their 
equations of motion at the same order in the form of ordinary differential equations fT^ (the 
result agrees with previous works Il^-|l9[); see Section 7. 



II. MULTIPOLE DECOMPOSITION 

In this section we construct the multipole expansion A4{h'^'^) = /i^^^t of the gravitational 
field outside an isolated system, supposed to be at once self-gravitating and slowly-moving. 
By slowly-moving we mean that the typical current and stress densities are small with respect 
to the energy density, in the sense that 
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where 1/c denotes (slightly abusively) the small post-Newtonian parameter. The point 
about (12) is that the ratio between the size of the source d and a typical wavelength of 
the gravitational radiation is of order d/X = 0(l/c). Thus the domain of validity of the 
post-Newtonian expansion covers the source: it is given by r < 6 where the radius b can be 
chosen so that d <b = 0{X/c). 



A. The matching equation 

The construction of the multipole expansion is based on several technical assumptions, 
the crucial one being that of the consistency of the asymptotic matching between the exterior 
and interior fields of the isolated system. In some cases the assumptions can be proved from 
the properties of the exterior field /ig^t as obtained in by means of a post-Minkowskian 



algorithm. However, since our assumptions are free of any reference to the post-Minkowskian 
expansion, we prefer to state them more generally, without invoking the existence of such an 
approximation (refer to p5| for the full detailed assumptions). In many cases the assump- 
tions have been explicitly verified at some low post-Newtonian orders [^,36,^,73 . 



The field h (skipping space-time indices) , solution in iR^ of the relaxed field equations and 
the no-incoming radiation condition, is given as the retarded integral (10). We now assume 
that outside the isolated system, say, in the region r > IZ where 7?. is a constant radius 
strictly larger than d, we have h = Ai{h) where Ai{h) denotes the multipole expansion of 
h, a solution of the vacuum field equations in IR^ deprived from the spatial origin r = 0, 
and admitting a spherical-harmonics expansion of a certain structure (see below). Thus, in 
Mx ml where Ml = - {0}, 

d^M{¥"') = 0, (2.2a) 
nMih^"") = MiA^") . (2.2b) 

The source term A^(A) is obtained from inserting Ai{h) in place of h into (5), i.e. Ai{A) = 
A{A4{h)). [Since the matter tensor has a compact support, A4(T) = so that A4{t) = 
Y^;^A^(A).] Of course, inside the source (when r < d), the true solution h differs from 
the vacuum solution Ai{h), the latter becoming in fact singular at the origin (r = 0). We 
assume that the spherical-harmonics expansion of Ai{h) in M x Ml reads 

M{h){^,t) = ^nir'^(lnr)Va,p(t) + i?7v(x,t) . (2.3) 

a<N 

This expression is valid for any G IN. The powers of r are positive or negative, a G 
and we have a < N (the negative powers of r show that the multipole expansion is singular 
at r = 0). For ease of notation we indicate only the summation over a, but there are two 
other summations involved: one over the powers p G IV of the logarithms, and one over 
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the order of multipolarity / G IN. The summations are considered only in the sense of 
formal series, as we do not control the mathematical nature of the series. The factor is 
a product of / unit vectors, rii = = n*i...n*', where L = is a multi-index with / 

indices, on which the symmetric and trace- free (STF) projection is applied: fii = STF[ni;,]. 
The decomposition in terms of STF tensors ni{6,ip) is equivalent to the decomposition in 
usual spherical harmonics. The functions LFa,p{t) are smooth {C°°) functions of time, which 
become constant when t < — T because of our assumption (9). [Of course, the L-^a,p's depend 
also on c : LFa,p{t, c).] Finally the function -RAr(x, t) is defined by continuity throughout IR^. 
Its two essential properties are Rj^ G C'^{IR^) and Rj^ = 0{r^) when r ^ with fixed t. In 
addition Rj^ is zero before the time — T. Though the function RN{-x,t) is given "globally" 
(as is the multipole expansion), it represents a small remainder 0{r^) in the expansion of 
Ai{h) when r — 0, which is to be identified with the "near-zone" expansion of the field 
outside the source. It is convenient to introduce a special notation for the formal near-zone 
expansion (valid to any order A^): 

Ap)(x,t) = 5^nir'^(lnr)^F„,p(t) , (2.4) 

where the summation is to be understood in the sense of formal series. [Note that (14) and 
(15) are written for the field variable A4{h), but it is easy to check that the same type of 
structure holds also for the source term A1(A).] 

Our justification of the assumed structure (14) is that it has been proved to hold for 
metrics in the class of Multipolar-post-Minkowskian (MPM) metrics considered in |23|, i.e. 
formal series /lext = ^"^n which satisfy the vacuum equations, are stationary in the past, 
and depend on a finite set of independent multipole moments. More precisely, from the 
theorem 4.1 in [^, the general MPM metric /lext, that we identify in this paper with A4{h'^'^), 
is such that the property (14) holds for the hnS to any order n, with the only difference that 
to any finite order n the integers a, p, I vary into some finite ranges, namely a^ninin) < a < N, 
< p < n — 1 and < Z < /max('T-), with amin(^) ~^ —oo and /max(^) — ^ +oo when n +oo. 
The functions L-F^p and the remainder i^^r in (14) should therefore be viewed as post- 
Minkowskian series ^G^iF^pn and J2^"'^N,n- What we have done in writing (14) and 
(15) is to assume that one can legitimately consider such formal post-Minkowskian series. 
Note that because the general MPM metric represents the most general solution of the field 



equations outside the source (Theorem 4.2 in |]23|), it is quite appropriate to identify the 



general multipole expansion Ai{h) with the MPM metric hext- Actually we shall justify this 
assumption in Section 5 by recovering from Ai{h), step by step in the post-Minkowskian 
expansion, the MPM metric /igxt- Because the properties are proved in [^ for any n, and 



because we consider the formal post-Minkowskian sum, we see that (14)-(15), viewed as if 
it were "exact", constitutes a quite natural assumption. In particular we have assumed in 
(14)- (15) that the multipolar series involves an infinite number of independent multipoles. In 
summary, we give to the properties (14)-(15) a scope larger than the one of MPM expansions 
(maybe they could be proved for exact solutions), at the price of counting them among our 
basic assumptions. 
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The multipole expansion A4{h) is a mathematical solution of the vacuum equations in 
M X Ml, but whose "multipole moments" (the functions LFa,p) are not determined in terms 
of the source parameters. When the isolated system is slowly moving in the sense of (12), 
there exists an overlapping region between the domains of validity of the post-Newtonian 
expansion: the "near- zone" r < b, where d < b = 0{X/c), and of the multipole expansion: 
the exterior zone r > TZ. For this to be true it suffices to choose TZ, which is restricted 
only to be strictly larger than d, such that d < IZ < b. We assume that the field h given 
by (10) admits in the near-zone a formal post-Newtonian expansion, h = h when r < b. 
On the other hand, recall that h = Ai{h) when r > TZ. Matching the two asymptotic 
expansions h and Ai{h) in the "matching" region TZ < r < b means that the (formal) 
double series obtained by considering the multipole expansion of (all the coefficients of) 
the post-Newtonian expansion h is identical to the double series obtained by taking the 
near-zone expansion of the multipole expansion. [We use the same overbar notation for the 
post-Newtonian and near- zone expansions because the near-zone expansion (r/c ^ 0) of the 
exterior multipolar field is mathematically equivalent to the expansion when c ^ oo with 
fixed multipole moments.] The resulting matching equation reads 



M{h)=M{h). (2.5) 

This equation should be true term by term, after both sides of the equation are re-arranged 
as series corresponding to the same expansion parameter. Though looking quite reasonable 
(if the theory makes sense), the matching equation cannot be justified presently with full 
generality; however up to 2PN order it was shown to determine a unique solution valid 



everywhere inside and outside the source |30,66|,|2^. The matching assumption complements 



the framework of MPM approximations by giving physical "pith" to the arbitrary 



multipole moments used in the construction of MPM metrics (see Section 4). 



B. The field in terms of multipole moments 

Let us consider the relaxed vacuum Einstein equation (13b), whose source term A^(A), 
according to our assumptions, owns the structure (14) [recall that (14) applies to J^{h) as 
well as A1(A)]. We obtain a particular solution of this equation (in M x Ml) as follows. First 
we multiply each term composing A^(A) in (14) by a factor (r/ro)^, where S is a complex 
number and ro a constant with the dimension of a length. For each term we can choose the 
real part of B large enough so that the term becomes regular when r — > 0, and then we can 
apply the retarded integral (11). The resulting 5-dependent retarded integral is known to 
be analytically continuable for any i? e C except at integer values including in general the 
value of interest B = 0. Furthermore one can show that the finite part (in short FPs=o) 
of this integral, defined to be the coefficient of the zeroth power of B in the expansion 
when i? ^ 0, is a retarded solution of the corresponding wave equation. In the case of a 
regular term in (14) such as the remainder Rn, this solution simply reduces to the retarded 
integral. Summing all these solutions, corresponding to all the separate terms in (14), we 
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thereby obtain as a particular solution of (13b) the object FP B=oOR^[ir /ro)^ J^{A)]. This 
is basically the method employed in to solve the vacuum field equations in the post- 
Minkowskian approximation. 

Now all the problem is to find the homogeneous solution to be added to the latter 
particular solution in order that the multipole expansion Ai{h) matches with the post- 
Newtonian expansion h, solution within the source of the field equation (3) [or, rather, 
(10)]. Finding this homogeneous solution means finding the general consequence of the 



matching equation (16). The result p4|,ill that the multipole expansion h'^'^ satisfying 



the Einstein equation (10) together with the matching equation (16) reads 

M{hn = FPb=o n-'[{r/rofM{An] - ^ E {I'^L^t - r/c)\ (2.6) 

where the first term is the previous particular solution, and where the second term is a re- 
tarded solution of the source-free (homogeneous) wave equation, whose "multipole moments" 
are given explicitly by {u = t — r/c) 

nl^iu) = FPb=o j d'^ |x/ro|^Xir'^^(x,M) . (2.7) 

Here T^'^ denotes the post-Newtonian expansion of the stress-energy pseudo-tensor r^'^ ap- 
pearing in the right side of (10). In (17) and (18) we denote L = ii . . . ii and Ol = di^ . . . 

Xi^ — Xi^ . . . Xij^ . 

It is important that the multipole moments (18) are found to depend on the post- 
Newtonian expansion T'^'^ of the pseudo-tensor, and not of t^^" itself, as this is precisely 
where our assumption of matching to the inner post-Newtonian field comes in. The formula 
is a priori valid only in the case of a slowly-moving source; it is a priori true only after in- 
sertion of a definite post-Newtonian expansion of the pseudo-tensor, where in particular all 
the retardations have been expanded when c —>■ oo [the formulas (17)-(18) assume implicitly 
that one can effectively construct such a post-Newtonian expansion]. 

Like in the first term of (17), the moments (18) are endowed with a finite part operation 
defined by complex analytic continuation in B. Notice however that the two finite part 
operations in the first term of (17) and in (18) act quite differently. In the first term of (17) 
the analytic continuation serves at regularizing the singularity of the multipole expansion at 
the spatial origin r = 0. Since the pseudo-tensor is smooth inside the source, there is no need 
in the moments (18) to regularize the field near the origin; still the finite part is essential 
because it applies to the bound of the integral at infinity (|x| — > oo). Otherwise the integral 
would be (a priori) divergent at infinity, because of the presence of the factor xl = 0{r^) in 
the integrand, and the fact that the pseudo-tensor t^^" is non-compact supported. The two 
finite parts present in the two separate terms of (17) involve the same arbitrary constant 
To, but this constant can be readily checked to cancel out between the two terms [i.e. the 
differentiation of M.{h^^) with respect to tq yields zero]. 

The formulas (17)-(18) were first obtained (in STF form) up to the 2PN order in |2^ by 



performing explicitly the matching. This showed in particular that the matching equation 
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(16) is correct to 2PN order. Then the proof vahd to any post-Newtonian order, but at the 
price of assuming (16) to all orders, was given in Section 3 of (see also Appendix A of 
for an alternative proof). The crucial step in the proof is to remark that the finite part 



of the integral of A4 (A) over the whole space is identically zero by analytic continuation: 

FPb=o [ rf'x |x/ro|^Xi>l(A)(x,u) = . (2.8) 



This follows from the fact that A^(A) can be written as a formal series of the type (15). Using 
(15) it is easy to reduce the computation of the integral (19) to that of the elementary radial 
integral J^°° d\x.\ |x|-^+2+'+'' (since the powers of the logarithm can be obtained by repeatedly 
differentiating with respect to B). The latter radial integral can be split into a "near-zone" 
integral, extending from zero to radius TZ, and a "far-zone" integral, extending from TZ to 
infinity (actually any finite non-zero radius fits instead of TZ). When the real part of S is a 
large enough positive number, the value of the near-zone integral is TZ^^^^'''^"' /{B + 3 + l + a), 
while when the real part of -B is a large negative number, the far-zone integral reads the 
opposite, —7^-^+^+'+"/ (5 + 3 + I + a). Both obtained values represent the unique analytic 
continuations of the near-zone and far- zone integrals for any i? G C except — 3 — / — a. The 
complete integral Jg^°° (i|x| |x|^+^+'+" is defined as the sum of the analytic continuations of 
the near- zone and far-zone integrals, and is therefore identically zero (Vi? G C); this proves 
(19). 

One may ask why the whole integration over contributes to the multipole moment 
(18) - a somewhat paradoxical fact because the integrand is in the form of a post-Newtonian 
expansion, and is thus expected to be physically valid (i.e. to give accurate results) only in 
the near zone. This fact is possible thanks to the technical identity (19) which enables us to 
transform a near- zone integration into a complete iR^-integration (refer to p5[ for details). 



C. Equivalence with the Will- Wiseman multipole expansion 



Recently a different expression of the multipole decomposition, with correlatively a differ- 
ent expression of the multipole moments, was obtained by Will and Wiseman |]7S[, extending 
previous work of Epstein and Wagoner and Thorne Basically, the multipole mo- 



ments in l?^ are defined by an integral extending over a ball of finite radius TZ (essentially 
the same TZ as here), and thus do not require any regularization of the bound at infinity. 
By contrast, our multipole moments (18) involve an integration over the whole M^, which 
is allowed thanks to the analytic continuation [leading to the identity (19)]. Let us outline 
the proof of the equivalence between the Will- Wiseman formalism and the present one 

\m- 

Will and Wiseman |7H] find, instead of (17)-(18), 



1=0 ' 



(2.9) 
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The first term is given by the retarded integral (11) acting on Ai{A), but truntated, as 
indicated by the subscript TZ, to extend only in the "far zone" : |x'| > 7^ in the notation (11). 
Thus, the near-zone part of the retarded integral, which contains the source, is removed, 
and there is no problem with the singularity of the multipole expansion at the origin. Then, 
the multipole moments Wl are given by an integral extending over the "near zone" only: 

Wr(^i)=/ d^XXLTf^'i^^u) . (2.10) 

The integral being compact-supported is well-defined. The multipole moments Wl look 
technically more simple than ours given by (18). On the other hand, practically speaking, 
the analytic continuation in (18) permits deriving many closed-form formulas to be used in 
applications [^,0 . Of course, one is free to choose any definition of the multipole moments 



as far as it is used in a consistent manner. 

We compute the difference between the moments 7{l and Wl- For the comparison we 
split T-Cl into far-zone and near- zone integrals corresponding to the radius TZ. Since the 
analytic continuation factor in T-Cl deals only with the bound at infinity, it can be removed 
from the near- zone integral, which is then clearly seen to agree with Wl- So the difference 
TCl — Wl is given by the far-zone integral: 

nM - Wl(m) = FPb=o / |x/ro| Vr(x, u) . (2.11) 

i|x|>7e 

Next we transform the integrand. Successively we write T = A4(t) because we are in the 
far zone; A4{t) = J^{t) from the matching equation (16); and Ai{T) = y^^A^(A) because 
T has a compact support. At this stage, the technical identity (19) allows one to transform 
the far-zone integration into a near zone integration (changing simply the overall sign in 
front of the integral). So, 

4 f 

HLiu) - Wl{u) = --^FPb=o / rf'x |x/ro| VMA)(x,^z) . (2.12) 

It is straightforward to check that the right side of this equation, when summed up over all 
multipolarities /, accounts exactly for the near-zone part that was removed from the retarded 
integral of A4{A) [first term in (20)], so that the "complete" retarded integral as given by 
the first term in (17) is exactly reconstituted. In conclusion the two formalisms [0,^] and 
75| are equivalent. 



III. SOURCE MULTIPOLE MOMENTS 

Quite naturally our source multipole moments will be closely related to the T^l's obtained 
in (18). However, before giving a precise definition, we need to find the equivalent of 
the multipole decomposition (17)-(18) in terms of symmetric and trace-free (STF) tensors, 
and we must reduce the number of independent tensors by imposing the harmonic gauge 
condition (13a). This leads to the definition of a "linearized" metric associated with the 
multipole expansion Ai{h), and parametrized by six sets of STF source multipole moments. 
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A. Multipole expansion in symmetric-trace-free form 



The moments TIl given by (18) are non-trace-free because xl owns all its traces (i.e. 
^iiii^x^L = where L — 2 = ii...ii_2). Instead of Ti^, there are certain advantages 

in using STF multipole moments: indeed the STF moments are uniquely defined, and they 
often yield simpler computations in practice. It is not difficult, using STF techniques, to 
obtain the multipole decomposition equivalent to (17)-(18) but expressed in terms of STF 
tensors. We find 

M{hn = FPs=oa-^'[{r/rorM{An] - ^ £ {I^'l' - r/c)\ (3.1) 

where the STF multipole moments are given by p5 | 

YVb=o j d?^\^/rQ\^XL j dz 5i{z)t>"'{-k,u + z\-x\/c) . (3.2) 



The notation for a STF product of vectors is = STF(a;2,) (such that is symmetric in 
L and Si^i^_^XL = 0; for instance Xij = XiXj — ^^jj-x^). As we see, the STF moments (25) 
involve an extra integration, over the variable z, with respect to the non-STF ones (18). 
The weighting function associated with the 2;-integration reads, for any I, 



(2/ + 1)!! 



l-z 



2\l . 



dz 6i{z) = 1 . 



(3.3) 



In the limit of large / the weighting function tends toward the Dirac delta measure (hence 
its name): lim;_»oo 5i = S. Remark that since (25) is valid only in the post-Newtonian 
approximation, the 2;-integration is to be expressed as a post-Newtonian series. Here is the 



relevant formula 30 



dz 6i{z)t{x, u + z|x| 



E 

k=0 



(2/+ 1) 



2'=A;!(2/ + 2A; + 1)! 



Ixl d , 



(3.4) 



In the limiting case of linearized gravity, one can neglect the first term in (24), and the 
pseudo-tensor t'^'^ in (25) can be replaced by the matter stress-energy tensor T^'^ (we have 
T^'^ = T^^ inside the slowly-moving source). Since T^^^ has a compact support the finite 
part prescription can be removed, and we recover the known multipole decomposition cor- 
responding to a compact-support source (see the appendix B in pO|). 



B. Linearized approximation to the exterior field 

Up to now we have solved the relaxed f\.e\d equation (10) in the exterior zone, with result 
the multipole decomposition (24)-(25). In this section we further impose the harmonic 
gauge condition (13a), and from this we find a solution of the linearized vacuum equation. 
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appearing as the first approximation in a post-Minkowskian expansion of tlie multipole 
expansion M.{h). 

Let us give a notation to tlie first term in (24): 



YVB=on-j,^[{r/r,)^M{kn] ■ (3.5) 



Applying on (24) the condition dyM.{h^^) = 0, we find that the divergence w'^ = d^u^^ is 
equal to a retarded solution of the source-free wave equation, given by 



'""^^(E^^^^^nt-r/M). (3.6) 



d' ^K^ l\ 

1=0 



Now, associated to any w'^ of this type, there exists some v'^'^ which is like w'^ a retarded 
solution of the source- free wave equation, n{v'^'') = 0, and furthermore whose divergence is 
the opposite of w^^, d^v'^" = —w^^. We refer to p3| , [70| for the explicit formulas allowing the 
"algorithmic" construction of v^"^ once we know w^. For definiteness, we adopt the formulas 
(2.12) in ||70[, which represent themselves a slight modification of the earlier formulas (4.13) 
in 1^ (see also the appendix B in [P5|]). 

With v'^'^ at our disposal we define what constitutes the linearized approximation to 
the exterior metric, say Gh^'^ where we factorize out G in front of the metric in order to 
emphasize its linear character: 

GhT = E ir^^ {l^L^^t - r/c)\ - v>^^ . (3.7) 

^ 1=0 ' ^ 

The linearized metric satisfies the linearized vacuum equations in harmonic gauge: Dh'^" = 
since both terms in (30) satisfy the source-free wave equation, and di,hi'^ = thanks to (29) 
and dyV^^ = —w^. Using the definition (30) one can re-write the multipole expansion of the 
exterior field as 

M {h^"") = Gh'^" + u>"' + v>"' . (3.8) 

Quite naturally the u^''' and v^'^ will represent the non-linear corrections to be added to the 
"linearized" metric G/i^"^ in order to reconstruct the complete exterior metric (see Section 
4)- ^ 

Since /ij"^ satisfies Dh'^'^ = = dyh^^ , there is a unique way to decompose it into the 
sum of a "canonical" metric introduced by Thorne |6^ (see also [^) plus a linearized gauge 
transformation, 

hT = KLi + d^vl + d^^'i - V'^'d.^l . (3.9) 
The canonical linearized metric is defined by 
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l>0 ' ^ 



1>1 ' ^ ^ ^ 



^ Y^eiahdaL-\(^JhL-\{u)^ , (3.10b) 

where the /^'s and J^'s are two sets of functions of the retarded time u = t — r/c [the 
subscript (n) indicates n time derivatives], and which are STF with respect to all their 
indices L = (the symmetrization is denoted with parenthesis). As for the gauge 

vector (f'^, it satisfies Uif'^ = and depends in a way similar to (33) on four other sets of 
STF functions of u, denoted Wl, Xl, Yl and Zl (one type of function for each component 
of the vector). See [ESI for the expression of ip^ = (p^Wl, Xl, Yl, Zl]. 



C. Derivation of the source multipole moments 

The two sets of multipole moments II and Jl parametrizing the metric (33) constitute 
our definitions for respectively the mass-type and current-type multipole moments of the 
source. Actually, there are also the moments Wl, Xl, Yl, Zl, and we refer collectively to 
{II, Jl, Wl, Xl, Yl, Zl} as the set of six source multipole moments. 

With (32) it is easily seen (because Dy?^ = 0) that the gauge condition d^h'^'^ = 
imposes no condition on the source moments except the conservation laws appropriate to 
the gravitational monopole I (having / = 0) and dipoles li, Ji {I = 1): namely, 

J(^) = 0; /f^ = 0; Ji'^ = 0. (3.11) 

The mass monopole / and current dipole Ji are thus constant, and agree respectively with 
the ADM mass and total angular momentum of the isolated system (later we shall denote 
the ADM mass by M = /). According to (34) the mass dipole Jj is a linear function of time, 
but since we assumed that the metric is stationary in the past, /, is in fact also constant, 
and equal to the (ADM) center of mass position. 

The expressions of II and Jl (as well as of the other moments Wl, Xl,Yl, Zl) come 
directly from (30) with (32)-(33) and the result of the matching, which is personified by the 
formula (25). To simplify the notation we define 

^00 _j_ jii 

S = ^ , (3.12a) 
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Si = — , (3.12b) 

c 

E,, = , (3.12c) 

(where = dijf^^). The result is 

/" R ( 4(2/ + 1) 

/l(m) = FPb=o y c?^x |x/ro| J dz^SiXL^, - ^2(-/ _^ ]^)(2/ + 3^ 

2(2/ + 1) ^ 
+ c^(/ + l)(/ + 2)(2/ + 5) ^'-^^^-^-"^^'^-^7^^^ """"/'^ ' ^^-^^"^ 



2/ + 1 



c2(/ + 2)(2/ + 3) 



^«+i^;L-i>ac5tS6c}(x,M + 2;|x|/c) , (3.13b) 



(<> refers to the STF projection). In a sense these expressions are exact, since they are 
formally valid up to any post-Newtonian order. [See (68)-(69) below for explicit formulas at 
2PN.] 

By replacing r^'^ in (36) by the compact-support matter tensor T'^'^ we recover the 
expressions of the multipole moments worked out in linearized gravity by Damour and Iyer 
|78| (see also |^). On the other hand the formulas (36) contain the results obtained by 



explicit implementation ( "order by order" ) of the matching up to the 2PN order 



IV. POST-MINKOWSKIAN APPROXIMATION 

In linearized gravity, the source multipole moments represent also the moments which are 
"measured" at infinity, using an array of detectors surrounding the source. However, in the 
non-linear theory, the gravitational source A^^" cannot be neglected and the first term in (24) 
plays a crucial role, notably it implies that the measured multipole moments at infinity differ 
from the source moments. Thus, we must now supplement the formulas of the source multi- 
pole moments (36) by the study of the "non- linear" term uf"" = FPB=QO^^[{r/ro)^M{A'"')] 



in (24). For this purpose we develop following ||2^ a post-Minkowskian approximation for 
the exterior vacuum metric. 



A. Multipolar-post-Minkowskian iteration of the exterior field 

The work started already with the formulas (31)-(33), where we expressed the exterior 
multipolar metric /ig^t = ■M.{h'^'^) as the sum of the "linearized" metric Gh^'^ and the "non- 
linear" corrections u^'^, given by (28), and v'^'^, algorithmically constructed from = d,^u^^ 
[see (29)]. The linearized metric is a functional of the source multipole moments: hi = 
hi[I, J, W, X, y, Z]. We regard G as the book-keeping parameter for the post-Minkowskian 
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series, and consider that Ghi is purely of first order in G, and thus that hi itself is purely 
of zeroth order. Of course we know from the previous section that this is untrue, because 
the source multipole moments depend on G; supposing hi = ©(G*^) is simply a convention 
allowing the systematic implementation of the post-Minkowskian iteration. 

Here we check that the non-linear corrections u'^'^ and v^'^ in (31) generate the whole post- 



Minkowskian algorithm of . The detail demanding attention is how the post-Minkowskian 
expansions of u^"^ and v'^'^ are related to a spliting of the gravitational source A'^'^ into 
successive non-linear terms. Let us pose, with obvious notation, 

A^" = N^'lh, h] + M^^'lh, h, h] + 0{h^) , (4.1) 

where, from the exact formula (5), the quadratic-order piece reads (all indices being lowered 
with the Minkowski metric, and h denoting rj^^hpcj): 

N^'^[h, h]= - hP''d%h^'' + ]^d^'hp,d''hP'' - ^d^^hd'h 

- d^^hp^d^h'^ - 9" V^^/i^'" + d^h^'Pid^h; + dph"^) 

--d^hp^d^hP" + \dphdPh + \dph,^d"h^ 

and where the cubic-order piece M[/i, /i, h] and all higher-order terms can be obtained in a 
straightforward way. 

First, reasoning ad absurdio, we prove (see for details) that both u and v indeed 



+ 1] 



(4.2) 



represent non-linear corrections to the linearized metric since they start at order G^: u = 
G^U2 + 0{G^) and v = G'^V2 + 0{G^). Next we obtain explicitly U2 by substituting the 
linearized metric hi into (38) and applying the finite part of the retarded integral, i.e. 

< = FFB=oa-n'{ir/ro)''N^''[hi,hi]} . (4.3) 

In this way we have a particular solution of the wave equation in iR x Ml, = N[hi, hi]. 
From U2 one deduces f 2 by the same "algorithmic" equations as used when deducing v from 
u [see after (29)]. Then □t>2 = and the sum U2 + V2 is divergenceless, so we can solve the 
quadratic-order vacuum equations in harmonic coordinates by posing 

/if = wf + v^" . (4.4) 

With this definition it is clear that the multipole expansion (31) reads to quadratic order: 

Mih^") = G/if + G'/if + 0{G'') . (4.5) 

Continuing in this fashion to the next order we find successively 

= FFB=oa-R'{r/rof(M^''[hi, hi, hi] + m^[hi, /ij + N^''[h2, /ii])}; 

(4.6a) 

/^r = < + ^3"^ (4.6b) 
Mih'^") = Gh^" + G^h^" + G'^h^" + 0{G^) . (4.6c) 
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This process continues ad infinitum. The latter post-Minkowskian algorithm is exactly the 
one proposed in (see also Section 2 of [0). That is, starting from hi[I, J,W,X,Y, Z] 
given by (32)- (33), one generates the infinite post-Minkowskian (MPM) series of solving 
the vacuum (harmonic-coordinate) Einstein equations in iR x Ml, and this formal series 
happens to be equal, term by term in G, to the general multipole decomposition of h^" 
given by (24). For any n, we have hl^'^ = u'^ + , and 

+ 00 

M{U"') = ^G''W^ . (4.7) 

n=l 

This result is perfectly consistent with the fact that the MPM algorithm generates the 
most general solution of the field equations in iR x M?^. Furthermore, the latter post- 



Minkowskian approximation is known [|26[ to be reliable (existence of a one-parameter family 
of exact solutions whose Taylor expansion when G ^ reproduces the approximation) - an 
interesting result which indicates that the multipole decomposition ]V[{K) given by (24)-(25) 
might be proved within a context of exact solutions. 

Recall that the source multipole moments /l, Jl, Wl-, ^l, Yl, Zl entering the linearized 
metric hi at the basis of the post-Minkowskian algorithm are given by formulas like (36). 
Thus, in the present formalism, the source moments, including formally all post-Newtonian 
corrections [and all possible powers of G] as contained in (36), serve as "seeds" for the 
post-Minkowskian iteration of the exterior field, which as it stands leads to all possible 
non-linear interactions between the moments. As we can imagine, rapidly the formalism 
becomes extremely complicated when going to higher and higher post-Minkowskian and/or 
post-Newtonian approximations. Most likely the complexity is not due to the formalism but 
reflects the complexity of the field equations. It is probably impossible to find a different 
formalism in which things would be much simpler (except if one restricts to a particular 
type of source). 



B. The "canonical" multipole moments 

The previous post-Minkowskian algorithm started with hi, a functional of six types of 
source multipole moments, II and Jl entering the "canonical" linearized metric /icani given 
by (33), and Wl, Xl, Yl, Zl parametrizing the gauge vector ipi in (32). All these moments 
deserve their name of source moments, but clearly the moments Wl, Xl, Yl and Zl do not 
play a physical role at the level of the linearized approximation, as they simply parametrize 
a linear gauge transformation. But because the theory is covariant with respect to (non- 
linear) diffeomorphisms and not merely to linear gauge transformations, these moments do 
contribute to physical quantities at the non-linear level. 

In practice, the presence of the moments Wl, Xl, Yl, Zl complicates the post- 
Minkowskian iteration. Fortunately one can take advantage of the fact (proved in [^]) that 
it is always possible to parametrize the vacuum metric by means of two and only two types 
of multipole moments Ml and Sl (different from II and Jl)- The metric is then obtained 
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by the same post-Minkowskian algorithm as in (39)-(43), but starting with the "canonical" 
linearized metric hca.ni[M, S] instead of hi[I, J, W, X, Y, Z]. The resulting non-linear metric 
/lean is isometric to our exterior metric h^xt = M.{h), provided that the moments Mi and 
Sl are given in terms of the source moments II, Jl, . . . , by some specific relations 

Ml = Ml[I,J,W,X,Y,Z] , (4.8a) 
Sl = Sl[I,J,W,X,Y,Z] . (4.8b) 

The two coordinate systems in which /lean and hext are defined satisfy the harmonic gauge 
condition in the exterior zone, but (probably) only the one associated with h^xt meshes with 
the harmonic coordinates in the interior zone. With the notation (32) the coordinate change 
reads Sx'^ = Gipi + non-linear corrections. We shall refer to the moments Ml and Sl as the 
mass-type and current-type canonical multipole moments. Of course, since at the linearized 
approximation the only "physical" moments are II and Jl, we have 

Ml = Il + 0{G), (4.9a) 
Sl = Jl + 0{G), (4.9b) 



where 0{G) denotes the post-Minkowskian corrections. Furthermore, it can be shown |73 
that in terms of a post-Newtonian expansion the difference between both types of moments 
is very small: 2.5PN order, i.e. 

ML = h + o(^^^ (4.10) 

[note that M = Madm = I]- Thus, from (46), the canonical moments are only "slightly" 
different from the source moments. Their usefulness is merely practical - in general they 
are used in place of the source moments to simplify a computation. 



C. Retarded integral of a multipolar extended source 

The previous post-Minkowskian algorithm has only theoretical interest unless we supply 
it with some explicit formulas for the computation of the coefficients hn- Happily for us 
pragmatists, such formulas exist, and can be found in a rather elegant way thanks to the 
process of analytic continuation. Basically we need the retarded integral of an extended (non- 
compact-support) source with a definite multipolarity /. Here we present three exemplifying 
formulas; see the appendices A in |]7D[ and [|7T| for more discussion. 

Very often we meet a wave equation whose source term is of the type nLF{t — r /c)/r^, 
where hi has multipolarity / and F denotes a certain product of multipole moments. [Clearly, 
the near- zone expansion of such a term is of the form (15).] When the power k is such that 
3 < A; < / + 2 (this excludes the scalar case / = 0), we obtain the solution of the wave 
equation as |p3|j68|l 
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[r/ro)^^F(t-r/c) 



X 



fc-3 
j=0 



k-3-j 



_ (fc-3)!(/ + 2-fc)! 
(/ + A;-2)! 

F^''-^-^\t-r/c) 



(4.11) 



As we see the (finite part of tlie) retarded integral depends in tliis case on tlie values of tlie 
extended source at tlie same retarded time t — rjc (for simplicity we use the same notation 
for the source and field points). But it is well known (see e.g. |]80|j8l| ) that this feature is 
exceptional; in most cases the retarded integral depends on the whole integrated past of 
the source. A chief example of such a "hereditary" character is the case with A; = 2 in the 
previous example, for which we find p8|,|69| 
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^F(t - r/c) 



ct~r 



/ dsF{s/c)Qi 



ct 



(4.12) 



where Qi denotes the Legendre function of the second kind, related to the usual Legendre 
polynomial Pi by the formula 



(4.13) 



Since the retarded integral (48) is in fact convergent when r — >■ 0, we have removed the factor 
{t/tq)^ and finite part prescription. When the source term itself is given by a "hereditary" 
expression such as the right side of (48), we get a more complicated but still manageable 
formula, for instance [TfT 
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(4.14) 



where denotes that ant i- derivative of F which is zero in the past [from (9) we have 

restricted F to be zero in the past], and where 



Rip{x) 



Q,{x) jyyQ,{y)'^{y) + Pi{x) jy dyQ,{y)'^{y) . (4.15) 



Like in (48) we do not need a finite part operation. The function Rip is well-defined thanks 
to the behaviour of the Legendre function at infinity: Qi{x) ~ l/x'"'^^ when x —>■ oo. 

The formulas (48)-(51) are needed to investigate the so-called tails of gravitational waves 
appearing at quadratic non-linear order, and even the tails generated by the tails themselves 
( "tails of tails" ) which arise at cubic order p^|7T| . (These formulas do not show a dependence 
on the constant tq, but other formulas do.) 



22 



V. RADIATIVE MULTIPOLE MOMENTS 



In Section 2 we introduced the definition of a set of multipole moments 
{II, Jl, Wl, Xl, Yl, Zl} for the isolated source, and in Section 3 we showed that the exte- 
rior field, and in particular the asymptotic field therein, is actually a complicated non-linear 
functional of the latter moments. Therefore, to define some source multipole moments is not 
sufficient by itself; this must be completed by a study of the relation between the adopted 
definition and some convenient far-field observables. The same is true of other definitions of 
source moments in different formalisms, such as in the Dixon local description of extended 
bodies []82|-^, which should be completed by a connection to the far- zone gravitational 
field, for instance along the line proposed by |85|j86[] in the case of the Dixon moments. In 
the present formalism, the connection rests on the relation between the so-called radiative 
multipole moments, denoted Ul and Vl, and the source moments II, Jl,- ■ ■, Zl [in fact, for 
simplicity's sake, we prefer using the two moments Ml and Sl instead of the more basic six 
source moments]. 



A. Definition and general structure 

The radiative moments Ul (mass-type) and Vl (current-type) are the coefficients of the 
multipolar decomposition of the leading part of the transverse-tracefree (TT) projection 
of the radiation field in radiative coordinates (T, X) (with i? = |X| the radial distance to 
the source). Radiative coordinates are such that the metric coefficients admit an expansion 
when i? — i> cx) in powers of 1/R (no logarithms of R). In radiative coordinates the retarded 
time T — R/c is light-like, or becomes asymptotically light-like when R —>■ oo. By definition, 



l>2 
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C{1 + Y)^'^L-2ecd(aVh)dL-2^ + O — 



(5.1) 



where iVj = R, Nl-2 = . . ■Ni^_^, Ncl-2 = NcNl-2, and the TT algebraic projector 
reads Vijab = (Sia - NiNa){5jb - NjNb) - ^^Sij - NiNj){6ab - NaNb). The radiative moments 
Ul and Vl depend on T — R/c; from (52) they are defined V / > 2. The radiative-coordinate 
retarded time differs from the corresponding harmonic-coordinate time by the well-known 
logarithmic deviation of light cones, 

^ R ^ r 2GM, fr\ . ^. 

T = t —ln{ — ]+0{G^), (5.2) 



c c \ro 



where we have introduced in the logarithm the same constant tq as in (39) (this corresponds 
simply to a choice of the origin of time in the far zone). 

Now from the post-Minkowskian algorithm of Section 3, it is clear that the radiative 
moments Ul and Vl can be obtained to any post-Minkowskian order in principle, in the 
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form of a non-linear series in the source or equivalently the canonical multipole moments 
Ml and Sl- The practical detail (worked out in |^9[) is to determine the transformation 
between harmonic and radiative coordinates, generalizing (53) to any post-Minkowskian 
order. The structure of e.g. the mass-type radiative moment is 

Ul = M« + 5] -^^;^,X^L . (5.3) 

n=2 

The first term comes from the fact that the radiative moment reduces at the linearized 
approximation to the (/th time derivative of the) source or canonical moment. The second 
term represents the series of non-linear corrections, each of them is given by a certain X^l 
which is a ra-linear functional of derivatives of multipole moments Ml or Sl- Furthermore 
we know from e.g. (48) and (50) that each new non-linear iteration (which always involves a 
retarded integral) brings a priori a new "hereditary" integration with respect to the previous 
approximation. So we expect that XnL is of the form (U = T — R/c) 

rU rU 

X„l{U)=J2 dm... du„Z,,{U,u,,...,Un)Mi^^'\u,)...Sil"\u„) (5.4) 

J ~oo J —oo 

where Zn denotes a certain kernel depending on time variables U,ui, . . . ,Un, and where 
the sum refers to all possibilities of coupling together the n moments. [See (56) below for 
examples of kernels Z2 and Z^.] A useful information is obtained from imposing that Zn be 
dimensionless; this yields the powers of G and 1/c in front of each non-linear term in (54), 
where k is the number of contractions among the indices present on the n moments (the 
current moments carrying their associated Levi-Civita symbol). 

As an example of application of (54) let us suppose that one is interested in the 3PN or 
1/c^ approximation. From (54) we have 3{n — 1) + 2k = 6, and we deduce that the only 
possibility is n = 3 (cubic non-linearity) and A; = (no contractions between the moments). 
From this we infer immediately that the only possible multipole interaction at that order is 
between two mass monopoles and a multipole, i.e. M x M x Ml- This corresponds to the 
"tails of tails" computed exphcitly in (56) below. 



B. The radiative quadrupole moment to 3PN order 

To implement the formula (54) a tedious computation is to be done, following in details 
the post-Minkowskian algorithm of Section 4 augmented by explicit formulas such as (47)- 
(51), and changing the coordinates from harmonic to radiative according to the prescription 
in . Here we present the result of the computation of the mass- type radiative quadrupole 
(/ = 2) up to the 3PN order: 

Um = MjfiU) + 2 V / drMjfiU - r) 



\2roJ 12 
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+- 



G 



dr 



a<i ]>a 



'J a<i ]>a 



2 /•+00 
'0 



+2(^)/ *M.f(C/-r) 



' 2ro/ 70 12^0/ 44100 



(5.5) 



Recall that in this formula the moment Mjj is the canonical moment which agrees with the 
source moment lij up to a 2.5PN term [see (46)], and that the source moment lij itself is 
given in terms of the pseudo-tensor of the source by (36a). See also the formulas (68)-(69) 
below for a more explicit expression of the source moment at the 2PN order [of course, to 
be consistent, one should use (56) conjointly with 3PN expressions of the source moments]. 

The "Newtonian" term in (56) corresponds to the quadrupole formalism. Next, there 
is a quadratic non-linear correction with multipole interaction M x M^- representing the 
dominant effect of tails (scattering of linear waves off the space-time curvature generated 
by the mass M). This correction, computed in is of order 1/c^ or 1.5PN and has the 
form of a hereditary integral with logarithmic kernel. The constant 11/12 depends on the 
coordinate system chosen to cover the source, here the harmonic coordinates; for instance the 
constant would be 17/12 in Schwarzschild-like coordinates |^^^. The next correction, of 
order 1/c^ or 2.5PN, is constituted by quadratic interactions between two mass quadrupoles, 
and between a mass quadrupole and a constant current dipole [Q. This term contains a 



hereditary integral, of a type different from the tail integral, which is due to the gravitational 
radiation generated by the stress-energy distribution of linear waves [p^- pll ; |69i ■ Sometimes 
this integral is referred to as the non-linear memory integral because it corresponds to the 



contribution of gravitons in the so-called linear memory effect [0. The non-linear memory 
integral can easily be found by using the effective stress-energy tensor of gravitational waves 
in place of the right side of (3); it follows also from rigorous studies of the field at future 
null infinity Finally, at 3PN order in (56) appears the dominant cubic non-linear 

correction, corresponding to the interaction M x M x M^- and associated with the tails of 



tails of gravitational waves [71 . 



C. Tail contributions in the total energy flux 

Observable quantities at infinity are expressible in terms of the radiative mass and current 
multipole moments. For instance the total gravitational-wave power emitted in all spatial 
directions (total gravitational flux or "luminosity" £) is given by the positive-definite mul- 
tipolar series 

^c2'+iU(/-l)/!(2/ + l)!!^ 
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4/(/ + 2) 



c2(/ _!)(/ + 1)1(2/ + 1) 



(5.6) 



In the case of inspiralling compact binaries (a most prominent source of gravitational waves) 
the rate of inspiral is fixed by the flux £, which is therefore a crucial quantity to predict. 
Excitingly enough, we know that C should be predicted to 3PN order for detection and 
analysis of inspiralling binaries in future experiments [p5|,p^. 

To 3PN order we can use the relation (56) giving the 3PN radiative quadrupole moment. 
Here we concentrate our attention on tails and tails of tails. The dominant tail contribution 
at 1.5PN order yields correspondingly a contribution in the total flux (with U = T — R/c): 



tail 



l\fiU) 



dTl\f{U-r) 



2rJ 12 



(5.7) 



Since we are interested in the dominant tail we have replaced using (46) the canonical 
mass quadrupole by the source quadrupole. Similarly there are some tail contributions due 
to the mass octupole, current quadrupole and all higher-order multipoles, but these are 
correlatively of higher post-Newtonian order [see the factors 1/c in (57)]. It has been shown 
||68| that the work done by the dominant "hereditary" contribution in the radiation reaction 
force within the source - which arises at 4PN order in the equations of motion - agrees 
exactly with (58). 

Next, because C is made of squares of (derivatives of) radiative moments, it contains 
a term with the square of the tail integral at 1.5PN. This term arises at the 3PN relative 
order and reads 



{tail)2 



5c 



11 



dTll'\U-T) 



In 



11 
12 



(5.8) 



Finally, there is also the direct 3PN contribution of tails of tails in (56): 
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(5.9) 



By a control of all the hereditary integrals in C up to 3PN we have checked that the 
terms (59)-(60) do exist. The two contributions (59) and (60) appear somewhat on the same 
footing - of course both should be taken into account in practical computations. Note that 
in a physical situation where the emission of radiation stops after a certain date, in the sense 
that the source multipole moments become constant after this date (assuming a consistent 
matter model which would do this at a given post-Newtonian order), the only contribution 
to C which survives after the end of emission is the 3PN tail-square contribution (59). 
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VI. POST-NEWTONIAN APPROXIMATION 



In Sections 2 and 3 we have reasoned upon the formal post-Newtonian expansion h^'^ 
of the near-zone field to obtain the source multipole moments as functionals of the post- 
Newtonian pseudo-tensor t^^'^. We have also considered in Sections 4 and 5 the formal 
expansion c oo of the radiation field when holding the multipole moments fixed. Clearly 
missing in this scheme is an explicit algorithm for the computation of h^'^ in the near zone. No 
such algorithm (say, in the spirit of the post-Minkowskian algorithm in Section 4) is known 
presently, but a lot is known on the first few post-Newtonian iterations [pO| , pl| ,^^- |58| , |30| , ^ 



The main difficulty in setting up a post-Newtonian algorithm is the appearance at some 
post-Newtonian order of divergent Poisson-like integrals. This comes from the fact that 
the post-Newtonian expansion is actually a near- zone expansion |^^, which is valid only in 
the region where r = 0(A/c), and that such an expansion blows up when taking formally 
the limit r —>■ +oo. For instance, Rendall has shown that the post-Newtonian solution 
cannot be asymptotically flat starting at the 2PN or 3PN level, depending on the gauge. 
This is clear from the structure of the exterior near-zone expansion (15), which involves many 
positive powers of the radial distance r. Thus, one is not allowed in general to consider the 
limit r +oo. In consequence, using the Poisson integral for solving a Poisson equation 
with non-compact-support source at a given post-Newtonian order is a priori meaningless. 
Indeed the Poisson integral not only extends over the near-zone but also over the regions at 
infinity. This means that the Poisson integral does not constitute the correct solution of the 
Poisson equation in this context. However, to the lowest post-Newtonian orders it works; for 
instance it was shown by Kerlick pO| , pT| and Caporali that the post-Newtonian iteration 
(including the suggestion by Ehlers of an improvement with respect to previous work 

is well-defined up to the 2.5PN order where radiation reaction terms appear, but that 
some divergent integrals show up at the 3PN order. 

Another difficulty is that the post-Newtonian approximation is in a sense not self- 
supporting, because it necessitates information coming from outside its own domain of 
validity. Of course we have in mind the boundary conditions at infinity which determine 
the radiation reaction in the source's local equations of motion. Again, to the lowest post- 
Newtonian orders one can circumvent this difficulty by considering retarded integrals that 
are formally expanded when c oo as series of "instantaneous" Poisson-like integrals |^ . 
However, this procedure becomes incorrect at the 4PN order, not to mention the problem 
of divergencies, because the near-zone field (as well as the source's dynamics) ceases to be 
given by an instantaneous functional of the source parameters, due to the appearance of 
"tail-transported" hereditary integrals modifying the lowest-order radiation reaction damp- 
ing 

Let us advocate here that the cure of the latter difficulty (and perhaps of all difficulties) 
is the matching equation (16). Indeed suppose that one knows a particular solution of the 
Poisson equation at some post-Newtonian order. This solution might be in the form of some 
"finite part" of a Poisson integral. The correct post-Newtonian solution will be the sum 
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of this particular solution and of a homogeneous solution satisfying the Laplace equation, 
namely a harmonic solution, regular at the origin, which can always be written in the form 
Y^A^xl, for some unknown constant tensors A^. The homogeneous solution is associated 
with radiation reaction effects. Now the matching equation states that the multipole expan- 
sion of the post-Newtonian solution agrees with the near- zone expansion of the exterior field 
(which has been computed beforehand in Section 4). The multipole expansion of the known 
particular solution can be obtained by a standard method, and the multipole expansion of 
the homogeneous solution is simply itself, i.e. M.{^AlXl) = YI^l^l- Therefore, we see 
that the matching equation determines in principle the homogeneous solution (i.e. all the 
unknown tensors Al), and since the exterior field satisfies relevant boundary conditions at 
infinity, the Al?> should correspond to the radiation reaction on a truly isolated system. See 
p7| , |68|j31| , |32[] for implementation of this method to determine the radiation reaction force to 
4PN order (1.5PN relative order). 



A. The inner metric to 2.5PN order 

Going to high post-Newtonian orders can become prohibitive because of the rapid prolif- 
eration of terms. Typically any allowed term (compatible dimension, correct index structure) 
does appear with a definite non-zero coefficient in front. However, high post-Newtonian or- 
ders can be manageable if one chooses some appropriate matter variables, and if one avoids 
expanding systematically the retardations due to the speed of propagation of gravity. Often 
it is sufficient, and clearer, to present a result in terms of matter variables still containing 
some c's, and perhaps also in terms of some convenient retarded potentials (being clear that 
any retardation going to an order higher than the prescribed post-Newtonian order of the 
calculation is irrelevant). See for instance (65) and (68)-(69) below. Anyway, only in a final 
stage, when a result to the prescribed order is in hands, should we introduce the more basic 
matter variables (e.g. the coordinate mass density) and perform all necessary retardations. 
Then of course one does not escape to a profusion of terms, but at least we have been able 
to carry the post-Newtonian iteration using some reasonably simple expressions. 



The matter variables are chosen pO| , p^ in a way consistent with our earlier definitions 
(35), i.e. 

CT = 5 ; (6.1a) 

(Ti = — ; (6.1b) 
c 

aij=r^ . (6.1c) 

To 2.5PN order one defines some retarded potentials V , Vi, Wij, X and Ri, with V and Vi 
looking like some retarded versions of the Newtonian and gravitomagnetic potentials, and 
Wij being associated with the matter and gravitational-field stresses: 
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V = D^^ {-inGa} , 



W, 



Ri 



□^^ <; -AirGiVcn, - Vicr) - 2dkVdiVk - -dtVdiV 



-AnGVau + 2VidtdiV + VdfV 



(6.2a) 
(6.2b) 
(6.2c) 

(6.2d) 



(6.2e) 



where denotes the retarded integral (11). All these potentials but V and Vi have a 
spatially non-compact support. The highest non-linearity entering them is cubic; it appears 
in the last term of X. 



Based on the latter potentials one can show |2^J7^ that the inner metric to order 2.5PN 



(in harmonic coordinates, 
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1 + -v^ + -v^' 1 + ^Wi, 
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(6.3a) 
(6.3b) 
(6.3c) 



(writing "g^^ would be more consistent with the notation of Section 2). With this form, we 
believe, the computational problems encountered in applications are conveniently divided 
into the specific problems associated with the computation of the various potentials (62), 
which constitute in this approach some appropriate computational "blocks" (having of course 
no physical signification separately). By expanding all powers of 1/c present into the matter 
densities (61) and into the retardations of the potentials (62), we find that the metric (63) 



becomes extremely complicated, as it really is (see e.g. P6| , ^ , |50| , |51| ) . 

Because of our use of retarded potentials, the metric (63) involves explicitly only even 
post-Newtonian terms (using the post-Newtonian terminology that even terms correspond 
to even powers of 1/c in the equations of motion). We have checked that the odd post- 
Newtonian terms (responsible for radiation reaction), contained in (63) via the expansion of 
retardations, match, in the sense of the equation (16), to the exterior metric satisfying the 
no- incoming radiation condition (9). 

The harmonic gauge condition implies some differential equations to be satisfied by the 
previous potentials. To 2.5PN order we find 



dAV + -, 



+ 2V^ 



Ri + VV, 



01^ 



(6.4a) 
(6.4b) 
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where Wa = 5ijWij. These equations are in turn equivalent to the equation of continuity 
and the equation of motion for the matter system, 




Note that the precision is IPN for the equation of continuity but only Newtonian for the 
equation of motion. 



From the 2.5PN metric (63) we obtain the pseudo-tensor r and the auxiliary quantities 
(35), that we replace into the formulas (36) to obtain the 2.5PN source multipole moments. 
Recall that the 2;-integration in the moments is to be carried out using the formula (27). 
Let us first see how this works at the IPN order. 

We need S to IPN order and Sj to Newtonian order. The latter quantity reduces to the 
matter part, Ej = + 0(l/c^), and the former one reads after a simple transformation 



The integrand is non-compact-supported because of the contribution of the second term, 
and accordingly we keep the regularization factor Ix/rol"^ and finite part operation. But let 
us operate by parts the second term, using the fact that |x|^XLA(y^) — A(|x|-^XL)y^ = 
di{\x\^XLdi{V'^) — di{\x\^XL)V'^} is a pure divergence. When the real part of S is a large 
negative number, we see thanks to the Gauss theorem that the latter divergence will not 
contribute to the moment, therefore by the unicity of the analytic continuation it will al- 
ways yield zero contribution. Thus, using Ax^ = 0, we can replace |x|-^a;iA(K^) in the 
second term of (67) by A{\:x.\^ x l)V'^ = B{B + / + 1)\x.\^^'^XlV'^ , and because of the explicit 
factor B we see that the second term can be non-zero only in the case where the factor B 
multiplies an integral owning a simple pole ~ l/S due to the integration bound |x| — > oo. 
Expressing (to Newtonian order) in terms of source points Zi and Z2, we obtain the 
integral J d^x |x|-^~^xl|x — Zi|~-'-|x — Z2\~^. When |x| — > oo each |x — Zi^2\~^ can be ex- 
panded as a series of ni,j j |x|~^i'2~^; then performing the angular integration shows that the 
sum of "multipolarities" I + li + I2 is necessarily an even integer. When this is realized 



B. The mass-type source moment to 2.5PN order 




(6.6) 



The substitution into the moments II given by (36a) leads to 




(6.7) 
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the remaining radial integral reads J (i|x| |x|^+'~'i~'2-2 ^j^ic]^ develops a pole only when 
I — li — I2 — 2 = —1. But that is incompatible with the previous finding. Thus the second 
term in (67) is identically zero, and we end up simply with a compact-support expression 
on which we no longer need to implement the finite part, 

4(2/ + l)xiL 



2c2(2/ + 3) ^* ^ ~ c2(/ + l)(2/ + 3)^*^^ ^ + ^ 



(6.8) 



This expression was first obtained in |^ using a different method valid at IPN order. Here 
we have recovered the same expression from the formula (36a) valid to any post-Newtonian 
order [^gg. 

Only starting at the 2PN order does the mass multipole moment have a non-compact 
support (so the finite part becomes crucial at this order). By a detailed computation in ^ 
we arrive at the following 2PN (or rather 2.5PN) expression: 
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(6.9) 



Recall that the canonical moment Ml differs from the source moment II at precisely the 
2.5PN order [see (46)]. 



VII. POINT-PARTICLES 

So far the post-Newtonian formalism has been developed for smooth (i.e. C°°) matter 
distributions. As such, the source multipole moments (36) become ill-defined in the presence 
of singularities. We now argue that the formalism is in fact also applicable to singular 
sources (notably point-particles described by Dirac measures) provided that we add to our 
other basic assumptions a certain method for removing the infinite self-field of point-masses. 
Our main motivation is the inspiralling compact binary - a system of two compact objects 
(neutron stars or black holes) which can be described with great precision by two point- 
particles moving on a circular orbit, and whose orbital phase evolution should be computed 
prior to gravitational-wave detection with relative 3PN precision P3|,PU . 



For this application we restrict ourselves to two point-masses nii and m2 (constant 
Schwarzschild masses). The trajectories are yi(t) and y2(t) and the coordinate velocities 
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vi,2 = dyi^2/dt] we pose = (c, vi^2)- The symbol 1^2 means the same term but 
with the labels of the two particles exchanged. A model for the stress-energy tensor of 
point-masses (say, at 2PN order) is 



T'"' fx f\ 

point— mass V^) 



^^,{t)v^,{t)vl{t)5[^-y,{t)] 



2; 



/ii(t) = 



mi 



(7.1a) 
(7.1b) 



v\ V 



where S denotes the three-dimensional Dirac measure, and g^^, the metric coefficients in 
harmonic coordinates {g = detg^,y). The notation [gg^u)! means the value at the location 
of particle 1. However, due to the presence of the Dirac measure at particles 1 and 2, the 
metric coefficients will be singular at 1 and 2. Therefore, we must supplement the model 
(70) by a method of "regularization" able to give a sense to the ill-defined limit at 1 or 
2. A priori the choice of one or another regularization constitutes a fully-qualified element 
of the model of point-particles. In the following we systematically employ the Hadamard 
regularization, based on the Hadamard "partie finie" of a divergent integral [p7| , p8| . 

Let us discuss an example. The "Newtonian" potential U, defined hj U = A^^(— 47rG'cr), 
where a is given by (61a) [we have V = U + 0(l/c^)], follows from (70a) as 
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(7.2) 



where ri = |x — yi|. To Newtonian order U = Gmi/ri + 0(1/0^) + 1 ^ 2. We compute U at 
the IPN order: from (70b) we deduce at this order /ii/mi = 1 — {U)i/c^ + vf/2c^ + 0{l/c^), 
which involves U itself taken at point 1, but of course this does not make sense because U 
is singular at 1 and 2. Now, after applying the Hadamard regularization (described below), 
we obtain unambiguously the standard Newtonian result {U)i = Gm2/ri2 + 0(l/c^), where 
^12 = |yi — y2|5 that we insert back into fii. So, U at IPN, and its regularized value at 1, 
read 
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A. Hadamard partie finie regularization 

We consider the class of functions of the field point x which are smooth on except at 
the location of the two source points yi_2, around which the functions admit some power-like 
expansions in the radial distance ri = |x — yi|, with fixed spatial direction ni = (x — yi)/ri 
(and idem for 2). Thus, for any -F(x) in this class, we have 
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^ = I]^i/i(a)(ni) (when n ^ 0) ; (7.4a) 

a 

F = Y1 ^2/2(a) M (when r2 ^ 0) , (7.4b) 

a 

where the summation index a ranges over values in Z bounded from below, a > — 
(we do not need to be more specific), and where the coefficients of the various powers of 
'"1,2 depend on the spatial directions rii 2. In (73) we do not write the remainders for the 
expansions because we don't need them; simply, we regard the expansions (73) as listings 
of the various coefficients /i(a) and /2(a)- We assume also that the functions F in this class 
decrease sufficiently rapidly when |x| 00, so that all integrals we consider are convergent 
at infinity. 

The integral J (P^F is in general divergent because of the singular behaviour of F near 



yi^25 but we can compute its partie finie (Pf) in the sense of Hadamard [^,0. Let us 
consider two volumes surrounding the two singularities, of the form ri < spi(ni) (and 
similarly for 2), where s measures the size of the volume and pi gives its shape as a function 
of the direction rii (pi = 1 in the case of a spherical ball). Using (73) it is easy to determine 
the expansion when s ^ of the integral extending on IB? deprived from the two previous 
volumes, and then to subtract from the integral all the divergent terms when s — in 
the latter expansion. The Hadamard partie finie is defined to be the limit when s of 
what remains. As it turns out, the result can be advantageously re-expressed in terms of 
an integral on IB? deprived from two spherical balls (pi 2 = 1), at the price of introducing 
two constants Si^2 which depend on the shape of the two regularizing volumes originally 
considered. With full generality the Hadamard partie finie of the divergent integral reads 

Pf / d^x F = limj / rf^x F 

J s^O{Jr,>s 

r2>s 

+ E / + 1^ (^) / ^^i/i(-3) + 1 2} (7.5) 
where si is given by 

Jdnj^.s) Inpi 

lnsi = — jr—— . (7.6) 

Because of the two arbitrary constants si,2 the Hadamard partie finie is ambiguous, and one 
could think a priori that there is no point about defining a divergent integral by means of 
an ambiguous expression. Actually the point is that we control the origin of these constants: 
they come from the coefficients of 2 in the expansions of F, which generate logarithmic 
terms in the integral. As we shall see the constants si^2 do not appear in the post-Newtonian 
metric up to the 2.5PN order (they are expected to appear only at 3PN order). 

We can also give a meaning to the value of the function F at the location of particle 
1 for instance, by taking the average over all directions rii of the coefficient of the zeroth 
power of ri in (73a), namely 
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(Fh - / ^/.(o, . (7.7) 

We refer also to the definition (76) as the Hadamard partie finie (of the function F at 1) 
because this definition is closely related to the definition (74) of the Hadamard partie finie 
of a divergent integral. To see this, apply (74) to the case where the function F is actually 
a gradient, F = diG, where G satisfies (73) [it is then clear that F itself satisfies (73)]. We 
find 

Pf J d^x diG = -A7iin\rlG)i - An{nirlG)2 (7.8) 

where in the right side the values at 1 and 2 are taken in the sense of the Hadamard partie 
finie (76). This nice connection between the Hadamard partie finie of a divergent integral 
and that of a singular function is clearly understood from applying the Gauss theorem on 
two surfaces ri 2 = s surrounding the singularities (there is no dependence on the constants 
■51,2)- 



B. Multipole moments of point-mass binaries 

To compute the source moments (36) of two point-particles we insert (70) in place of the 
stress-energy tensor T'^'^ of a continuous source, and we pick up the Hadamard partie finie 
[in the sense of (74)] of all integrals. This ansatz reads 

(-^L)point-mass = Pf {-^L[^point-mass] } i (7.9a) 
(</L)point-mass = Pf {^^[^^omt-mass] } • (7.9b) 

As we have seen in (69), the source multipole moments involve at high PN order many 
(non-compact-support) non-linear contributions which can be expressed in terms of retarded 
potentials such as V . The paradigm of such non-linear contributions is a term involving the 
quadratic product of two (derivatives of) potentials V , say dVdV, or, neglecting 0(l/c^) 
corrections, dUdU. To Newtonian order U is given by Grrii/ri + Gm2/r2 and it is easily 
checked that this paradigmatic term can be written as a certain derivative operator, say dd, 
acting on the elementary integral (assuming for simplicity / = 2) 

Y^,{yl,y2) ^ -^FPb=o [ ^'x |x/ro|^^ . (7.10) 

We see that the integral would be divergent at infinity without the finite part operation. 
However, it is perfectly well-behaved near 1 and 2 where there is no need of a regularization. 
The integral (79) can be evaluated in various ways; the net result is P^JT^ 
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(7.11) 
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where < ij >= STF(zj). Starting at 3PN order we meet some elementary integrals which 
need the regularization at 1 or 2 in addition to involving the finite part at infinity. An 
example is 



%(yi) = -i-Pf|FPB=oy" rf'x |x/ro|«^} . 



(7.12) 



To obtain this integral one splits it into a near-zone integral extending over the domain 
ri < TZi (say), and a far- zone integral extending over TZi < ri. The Hadamard regularization 
at 1 applies only to the near-zone integral, while the finite part at -B = is needed only for 
the far- zone integral. The result, found to be independent of the radius TZi, reads [|77| 

16- 
' 15- 

In this case we find an explicit dependence on both the constants tq due to the finite part at 
infinity, and Si due to the Hadamard partie finie near 1 [see (74)]. However these constants 
do not enter the multipole moments before the 3PN order (collaboration with Iyer and 
Joguet |]77|). 
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(7.13) 



A long computation, done in ||72[, yields the mass- type quadrupole moment at the 2PN 
order fully reduced in the case of two point-masses moving on a circular orbit. The method 
is to start from (69) (issued from |^) and to employ notably the elementary integral (79)- 
(80) (see also [|72[ for the treatment of a cubically non-linear term). An equivalent result has 
been obtained by Will and Wiseman using their formalism [^. In a mass-centered frame 
the moment is of the form 

■ 1 



(7.14) 



where i/i = y\—y\ and Vi = v{—V2, where u denotes the binary's Newtonian orbital frequency 
[u;^ = Gm/r\2 with m = mi +^2], and where fi = mim2/m is the reduced mass. The point 
is to obtain the coefficients A and B developed to 2PN order in terms of the post-Newtonian 
parameter 7 = Gm/r 120^, where we recall that ri2 is the distance between the two particles 
in harmonic coordinates. Untill 2PN we find some definite polynomials in the mass ratio 
1/ = fi/m (such that < z/ < 1/4): 
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(7.15b) 
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The 2PN mass quadrupole moment (83)- (84) is part of a program aiming at computing 
the orbital phase evolution of inspiralling compact binaries to high post-Newtonian order 
(see Section 7.4). First-order black- hole perturbations, valid in the test-mass limit i/ ^ 
for one body, have already achieved the very high 5.5PN order |]87|,|99|-[101[] . Recovering 



the result of black-hole perturbations in this limit constitutes an important check of the 
overall formalism. For the moment it passed the check to 2.5PN order [^,^; this is quite 
satisfactory regarding the many differences between the present approach and the black-hole 
perturbation method. 
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C. Equations of motion of compact binaries 



The equations of motion of two point-masses play a crucial role in accounting for the 
observed dynamics of the binary pulsar PSR1913-I-16 1)11-^,^^, and constitute an important 
part of the program concerning inspiralling compact binaries. The motivation for inves- 
tigating rigorously the equations of motion came in part from the salubrious criticizing 
remarks of Jiirgen Ehlers et al 0. Four different approaches have succeeded in obtaining 
the equations of motion of point-mass binaries complete up to the 2.5PN order (dominant 
order of radiation reaction): the "post-Minkowskian" approach of Damour, Deruelle and 
colleagues |]T6|-[ 



]; the "Hamiltonian" approach of Schafer and predecessors |p.02| , |103| . |57| . |58 
; the "extended-body" approach of Kopejkin et al ||104| , |105|] ; and the "post-Newtonian" ap- 
proach of Blanchet, Faye and Ponsot [|7^. The four approaches yield mutually agreeing 
results. 

The post-Newtonian approach consists of (i) inserting the point-mass stress-energy 
tensor (70) into the 2.5PN metric in harmonic coordinates given by (63); (ii) curing systemat- 
ically the self-field divergences of point-masses using the Hadamard regularization; and (iii) 
substituting the regularized metric into the standard geodesic equations. For convenience 
we write the geodesic equation of the particle 1 in the Newtonian-like form 



dt 



(7.16) 



where the (specific) linear momentum V\ and force T\ are given by 



r-- 



(7.17) 



Crucial in this method, the quantities are evaluated at the location of particle 1 according to 
the rule (76). All the potentials (62) and their gradients are evaluated in a way similar to our 
computation of U in (72), and then inserted into (85)-(86). We "order-reduce" the result, 
i.e. we replace each acceleration, consistently with the approximation, by its equivalent in 
terms of the positions and velocities as given by the (lower-order) equations of motion. After 
simplication we find, in agreement with other methods. 
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[where = {y{ — 1/2) /r^] = v\ — v^'-, and e.g. (r;,i2fi) denotes the Euchdean scalar 
product]. At the IPN or 1/c^ level the equations were obtained before by Lorentz an Droste 
and by Einstein, Infeld and Hoffmann [^. The 2.5PN or l/c' term represents the 



radiation damping in harmonic coordinates [correct because the metric (63) we started with 
matches to the post-Minkowskian exterior field] . In the case of circular orbits, the equations 
simplify drastically: 
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where the orbital frequency co'2PN of the 2PN circular motion reads 
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(the post-Newtonian parameter is 7 = Gm/c^r 12; and v = fi/m). 



D. Gravitational waveforms of inspiralling compact binaries 

The gravitational radiation field and associated energy flux are given by (52) and (57) 
in terms of time-derivatives of the radiative multipole moments, themselves related to the 
source multipole moments by formulas such as (56). Furthermore, at a given post-Newtonian 
order, the source moments admit some explicit though complicated expressions such as 
(68)-(69), which, when specialized to (non-spinning) point-mass circular binaries, yield e.g. 
(83)-(84). 

Now, for insertion into the radiation field and energy flux, one must compute the time- 
derivatives of the binary moments, with appropriate order-reduction using the binary's equa- 
tions of motion (87)- (89). This yields in particular the fully reduced (up to the prescribed 
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post-Newtonian order) gravitational waveform of the binary, or more precisely the two in- 
dependent "plus" and "cross" polarization states /i+ and /ix- The result to 2PN order is 
written in the form 

2Gmvx 



cm 



{Hfl + xV^i/jV^) + ^H^l + :,^nHf!:^ + x^H^l } , (7.21) 



where, for convenience, we have introduced a post-Newtonian parameter which is directly 
related to the orbital frequency: x = {Gmuj2PN / c^Y^^ , where u;2pn is given for circular orbits 
by (89). The various post-Newtonian coefficients in (90) depend on the cosine and sine of 
the "inclination" angle between the detector's direction and the normal to the orbital plane 
(cj = cost and Si = sini), and on the masses through the ratios u = fi/m and 6m/ m, where 
6m = mi — m2. The result for the "plus" polarization (collaboration with Iyer, Will and 
Wiseman 111 0611) is 
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The "cross" polarization admits a similar expression (see ||106|| ). Here, if) denotes a particular 
phase variable, related to the actual binary's orbital phase and frequency u = uj2pn by 
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is the angle, oriented in the sense of the motion, between the vector separation of the 
two bodies and a fixed direction in the orbital plane (since the bodies are not spinning, the 
orbital motion takes place in a plane). In (92), i^o denotes some constant frequency, for 
instance the orbital frequency when the signal enters the detector's frequency bandwidth; 
1U6II for discussion. 



see 



The previous formulas give the waveform of point-mass binaries whenever the frequency 
and phase of the orbital motion take the values uj and 0. To get the waveform as a function 
of time, we must replace uj and (j) by their explicit time evolutions uj{t) and (j){t). Actually, 
the frequency is the time-derivative of the phase: uj = d(f)/dt. The evolution of the phase is 
entirely determined, for circular orbits, by the energy balance equation dE/dt = —C relating 
the binding energy E of the binary in the center of mass to the emitted energy flux C. E 
is computed using the equations of motion (87), and £ follows from (57) and application of 
the previous formalism [changing the radiative moments to the source moments, applying 
(83)-(84), etc...]; the net result for the 2.5PN orbital phase [^,^|7| is 
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where 0o is a constant phase (determined for instance when the frequency is cuq), and G the 
convenient dimensionless time variable 

e = ^(te-t), (7.25) 

tc being the instant of coalescence at which, formally, uj{t) tends to infinity (of course, the 
post-Newtonian method breaks down before the final coalescence). All the results are in 
agreement, in the limit z/ — > 0, with those of black-hole perturbation theory P7|,p5|-[T(7T[ . 



VIII. CONCLUSION 

The formalism reviewed in this article permits investigating in principle all aspects of 
the problem of dynamics and gravitational- wave emission of a slowly-moving isolated system 
(with, say, u/c ~ 0.3 at most): the generation of waves, their propagation in vacuum, the 
back-reaction onto the system, the structure of the asymptotic field, and most importantly 
the relation between the far-field and the source parameters. Of course, the formalism is 
merely post-Newtonian and never "exact", but in applications to astrophysical objects such 
as inspiralling compact binaries this should be sufficient provided that the post-Newtonian 
approximation is carried to high order. 

Furthermore, there are several places in the formalism where some results are valid 
formally to any order of approximation. For instance, the source multipole moments are 
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related to the infinite formal post-Newtonian expansion of the pseudo-tensor [see (18) or 
(36)], and the post-Minkowskian iteration of the exterior field is performed to any non-linear 
order [see (43)]. In such a situation, where an infinite approximate series can be defined, 
there is the interesting question of its relation to a corresponding element in the exact theory. 
For the moment the only solid work concerns the post-Minkowskian approximation of the 
exterior vacuum field, which has been proved to be asymptotic [^. Likewise it is plausible 
that the expressions of the source multipole moments could be valid in the case of exact 
solutions. 

The most important part of the formalism where a general prescription for how to proceed 
at any approximate step is missing, is the post-Newtonian expansion for the field inside the 
isolated system. For instance, though the multipole moments are given in terms of the 
formal post-Newtonian expansion of the pseudo-tensor, no general algorithm for computing 
explicitly this post-Newtonian expansion is known. An interesting task would be to define 
such an algorithm, in a manner similar to the post-Minkowskian algorithm in Section 4. In 
the author's opinion, the post-Newtonian algorithm should be defined conjointly with the 
post-Minkowskian algorithm, and should rely on the matching equation (16), so as to convey 
into the post-Newtonian field the information about the exterior metric. 

Note that even if a general method for implementing a complete approximation series 
is defined, this method may be unworkable in practical calculations, because not explicit 
enough. For instance the post-Minkowskian series (43) is defined in terms of "iterated" 
retarded integrals, but needs to be suplemented by some formulas, to be used in applications, 
for the retarded integral of a multipolar extended source. In this respect it would be desirable 
to develop the formulas generalizing (50)-(51) to any non-linear order. This should permit 
in particular the study of the general structure of tails, tails of tails, and so on. 

For the moment the only application of the formalism concerns the radiation and motion 
of point-particle binaries. Of course it is important to keep the formalism as general as 
possible, and not to restrict oneself to a particular type of source, but this application to 
point-particles offers some interesting questions. Indeed, it seems that the post-Newtonian 
approximation used conjointly with a regularization d la Hadamard works well, and that one 
is getting closer and closer to an exact (numerical) solution corresponding to the dynamics 
and radiation of two black-holes. So, in which sense does the post-Newtonian solution 
(corresponding to point-masses without horizons) approach a true solution for black-holes? 
Does the adopted method of regularizing the self-field play a crucial role? Is it possible to 
define a regularization consistently with the post-Newtonian approximation to all orders? 
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